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The solar photospheric abundance of zirconium 
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Zirconium (Zr), together with strontium and yttrium, is an important element in the understanding of the Galactic nucle- 
osynthesis. In fact, the triad Sr-Y-Zr constitutes the first peak of s-process elements. Despite its general relevance not many 
studies of the solar abundance of Zr were conducted. We derive the zirconium abundance in the solar photosphere with the 
same C05B0LD hydrodynamical model of the solar atmosphere that we previously used to investigate the abundances 
of C-N-O. We review the zirconium lines available in the observed solar spectra and select a sample of lines to determine 
the zirconium abundance, considering lines of neutral and singly ionised zirconium. We apply different line profile fitting 
strategies for a reliable analysis of Zr lines that are blended by lines of other elements. The abundance obtained from lines 
of neutral zirconium is very uncertain because these lines are commonly blended and weak in the solar spectrum. How- 
ever, we believe that some lines of ionised zirconium are reliable abundance indicators. Restricting the set to Zrll lines, 
from the C05B0LD 3D model atmosphere we derive A(Zr)=2.62 ± 0.06, where the quoted error is the RMS line-to-line 
scatter. 
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1 Introduction 

The photospheric solar abundances of some elements have 
been studied more often than others; this is mainly due to 
the importance of the elements to explain nucleosynthesic 
processes. Moreover, the difficulty of extracting suitable 
lines with accurately known transition probabilities in the 
visible range of the solar spectrum can also explain why 
some elements have been studied less extensively. 

For the triad Sr-Y-Zr, there exist only a few detailed 
studies of their solar abundances (an ADSLj search yields 4 
papers for Sr, 5 for Y, and 10 for Zr) in spite of their impor- 
tance in Galactic chemical evolution. Sr-Y-Zr, Ba-La, and 
Pb are located at three abundance peaks of the s-processes 
producing the enrichment of these elements in the Galaxy. 
Knowledge of the present Sr-Y-Zr abundances in stars of 
different metallicity and age is required to understand the 
co mplicated nucleosynth esis of these elements (for details 
see 



Travagho et al]|2004 



Zirconium is present in the solar spectrum with lines 
of Zrl and Zrll. The dominant species is Zrll. Accord- 
ing to our 3D model, and in agreement with our ID ref- 
erence model, about 99% of zirconium is singly ionised 
in the solar photosphere. Departure from local thermo- 
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dynamical equilibrium (LTE) probably affects Zrl lines, 
producing a too low Zr abundance under the assumption 
of LTE. In fact NLTE effects, even on weak unsa turated 
Zrl lines, have been found by jBrownet al. ' ('19831) for G 



and K giants. However, both iBiemont et al . ( 1981,) and 
Bogdanovich et al.l ( 19961) analysed in LTE Zri lines (34 
and 21 lines, respectively) and Zrll lines (24 and 15 lines, 
respectively) in the solar photosphere, and found an excel- 
lent agreement between the abu ndances derived from both 
ionisation stages. Very recently Velichko et al. (2010l) per- 



formed a NLTE analysis of zirconium in the case of the 
Sun and late type stars. According to their computations 
the NLTE abundance is larger than the LTE one, by up to 
0.03 dex for Zr II and by 0.29 dex for Zr I. 

Owing to th e small number of Zr I and Zrll lines in 
metal poor stars ( Gratton & Snedenll994 ), it is important to 
analyse both of them in the solar photosphere to derive the 
solar abundance from both ionisation stages and to assess 
the agreement between the derived abundances. In spectra 
of cool stars it is eas i er to observe Zrl lines. For exam- 
ple, Goswami & Aoki ( 2010l) realised that none of the Zrll 
lines were usable in their analysis of the cool Pop. II CH star 
HD 209621, and the Zr abundance is derived from the only 
Zrl line in their spectrum at 61 3.457 nm. A similar situa- 
tion had been encountered bv lVanture & WallersteinI (120021) 
in their study of the Zr/Ti abundance ratio in cool S stars, 
where only Zr I lines in the red part of the spectrum can be 
used. 
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2 Lines and atomic data 



According to lMalcheva et al. I (l2006l) . zirconium has five sta- 
ble isotopes. Four of them (^°Zr, ^^Zr, ^^Zr, and ^''Zr) are 
produced by the s-process; the fifth, ^^Zr, is produced in the 
r-process. Zirconium is a very refractory element, it is dif- 
ficult to be vaporised by conventional thermal means and, 
consequently, has been relatively little investigated in the 
laboratory. 

In Table[T] and Table|2] the Zr I and Zr II lines, with the 
log gf values available in the literature are collected. 

iBiemont et al. ( 1981 ) selected lines that lie at A < 
80 nm. They derived log g f using the technique developed 
bv lHannaford & Lowel(ll981.) . suitable for highly refractory 
elements like Zr, to determine lifetimes of 34 levels of Zr I 
and 20 levels of Zrll. From these measurements, coupled 
with the measurements of branching ratios, they derived the 
log gf values of 38 Zr I and 3 1 Zr II lines. 

iBogdanovich et al" ( 1996 ) computed and used log gf 
values for 21 Zr I lines. They give ne w log gf values for 1 5 
Zr II lines among those selected by iBiemont et al. I (ll98lh . 
Their log gf are given in Table|2l 

Sikstrom et al. ( 19991) derived the abundance of zirco- 



nium in HgMn star x Lupi, finding a disagreement when 
using Zrll or Zrlll lines. They measured the f- values for 
several Zr II lines in the UV at A < 300 nm. 



Vanture & WallersteinI (120021) extended the search for 



Zr I fines in the near-IR, at wavelengths longer than 800 nm, 
a region important in the study of cool stars, because this 
is the region were cool stars emit most of their flux. The 
complete sample of these near-infrared lines is also given in 
Table[T] Because absorption bands of molecules are weak 
in the warmer S stars, these near-IR zirconium lines are bet- 
ter abundance indicators than the Zr II lines in the blue part 
of the spectrum. These lines are, however, not necessarily 
good for the Sun. 

We have al s o che cked t he 1 3 Zr i lines i dentifi ed by 
Swensson et"aD (Il970l) in the lOelbouiUe et al.1 (Il98lb atlas, 
but 12 of them appear severely blended. O nly the line at 



784. 9 nm i s retained in Tableffl as it was by IBiemont et al. 
(Il98lb and lVanture & WaUersteinl(l2002l) . 



In the present an alysis we adopt the log gf values of 
Biemont et all (Il98lb for all Zr I lines. 

According to the NIST database, Zrll lines lie only 
in the near UV-blue (241. 941-535.035 nm); only two weak 
lines (at 667.801 and 678.715 nm) are present in NIST with 
A >540nm. 

The line identifications by 'Moore et al. ' (Il966l) on the 
Utrecht solar atlas have been used by Biemo nt et al.l (Il98lh 
to select 31 Zr II lines, seven of which were later discarded 
because the derived abundance exceeded the mean by more 
than 3a , so they are likely blend ed with other unknown 
sp ecies. Gratton & Sneden ( 19941) used the same line list 



as 



Biemont et al.l (il981i) . ILiung et al.l (l2006h derived new 



oscillator strengths for 263 Zr II lines and studied 7 lines, 
the lines that they judged to be the best and unperturbed in 



the photospheric spectrum, to derive the solar Zr abundance 
based on both ID and 3D model atmospheres. 



We extracted from the 243 li nes by iMalcheya et al 



Biemont et al.l 



ilMll); 



(l2006h . those in common with 
for th ese lines the b est values of log gf, according to 
Malch eva et al. (2006), are the same as in IBiemont et al 



(1981). We compared the new log gf values determine d 
by iLjunget al.l(l2006l) with those by IBiemont et all (Il98ll) . 
There is a good agreement for most of the lines, with a few 
exceptions (see Table|2]i. We report all these values in Ta- 
blelH 

In the present analysis we used the log gf values of 
Ljung et al.l (l2006l) for all the Zr II lines. 



3 Solar Zr abundance in the literature 

Several analyse s of the solar Zr abunda nce, mad e be- 
fore 1980, used 'Corliss & Bozman ' (1962) log gf (' AUed 
1 1965 Wallersteini , 1 966, Grevesse et al. .1968.) or older ones 
dGoldberg et al.lll960l) or have been made to derive a tem- 
perature correction for the ICorliss & Bozman 1(11962') data 
tAUeni il976,) . The use of these transition probabilities, 
which are known to contain errors depending on excitation 
and temperature, coupled with the use of old solar spectral 
atlases, explains discrepancies of these Zr abundance deter- 
minations based either on Zr I and/or Zr II. 

We concentrate on the more recent determinations of the 
photospheric solar Zr abundance, that we summarise below. 



1. Biemon t et al.l (119811) analysed the iDelbouille et al 



(1973) atlas. The equivalent widths (EWs) were in- 
dependently measured by E. Biemont and N. Grevesse, 
and the results examined and discussed to produce the 
final list. T hey used the Holweger-M iiller solar model 
dHolwegerll 19671: iHolweger & MuUeiil 19741) wifli a mi- 
croturbulence, ^, of 0.8kms~^. They note that their 
analysis is independent of line-broadening parameters 
since most of the lines are faint or very faint. The 
log gf are from their own measurements. They find 
that the abundance from Zrl lines is independent of 
f and that the abundance from Zrll lines is insensi- 
tive to the model. These authors rejected 4 Zrl lines 
and 7 Zr II lines from their selected sample in their fi- 
nal solar abundance analysis because they imply a Zr 
abundance which is more than 3a higher than the aver- 
age. These lines appear after the horizontal line in Ta- 
bles[T] and |2] They conclude that Zr II lines are better 
abundance indicators. They obtain a zirconium abun- 
dance of A(Zr)=2.57 ± 0.07 from Zrl (34 fines) and 
A(Zr)=2.56 ± 0.07 from Zrll (24 lines). 
2. Gratton & SnedenI (Il99 4') studied the Zr behaviour in 
metal -poor stars . They used as model atmosphere a grid 



from iBell et al.l (119761) for giants, and a similar grid 
for dwarfs provided by Bell. They also analysed the 
solar spectrum, to have a reference A( Zr)0. For this 



purpose they used the same 24 lines as IBiemont et al 
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Table 1 Lines of Zr i chose n bv lBiemont et alj (119811) and 
Vanture & WallersteinI ( l2002h 



Ta ble 2 Lines of Z r ii chosen by iBiemont et alj (1198 ih 
and lLjung et al.l(l2006l) 



A 


Elow 


EW 




log gf 




nm 


eV 


pm 


B 


vw 


Bog 


350.9331 


0.07 


0.65 


-0.11 




-0.21 


360.1198 


0.15 


1.3 


-0.47 






389.1383 


0.15 


1.7 


-0.10 






402.893 


0.52 


0.06 


-0.72 






403.0049 


0.60 


0.26 


-0.36 




-0.59 


404.3609 


0.52 


0.58 


-0.37 






407.2696 


0.69 


0.57 


+0.31 




+0.24 


424.1706 


0.65 


0.37 


+0.14 




+0.07 


450.7100 


0.54 


0.36 


-0.43 




-0.46 


454.2234 


0.63 


0.46 


-0.31 






468.7805 


0.73 


1.00 


+0.55 




+0.30 


471.0077 


0.69 


1.05 


+0.37 




+0.19 


473.2323 


0.63 


0.25 


-0.49 




-0.56 


473.9454 


0.65 


0.55 


+0.23 




+0.07 


477.2310 


0.62 


0.53 


+0.04 




-0.07 


478.494 


0.69 


0.16 


-0.49 






480.589 


0.69 


0.15 


-0.42 




-0.63 


480.9477 


1.58 


0.16 


+0.16 






481.5056 


0.65 


0.20 


-0.53 




-0.22 


481.5637 


0.60 


0.30 


-0.03 






482.806 


0.62 


0.19 


-0.64 






504.655 


1.53 


0.050 


+0.06 




-0.25 


538.5128 


0.52 


0.18 


-0.71 






612.746 


0.15 


0.21 


-1.06 


-1.06 


-0.87 


613.457 


0.00 


0.19 


-1.28 


-1.28 


-1.05 


614.046 


0.52 


0.073 


-1.41 


-1.41 


-0.85 


614.3183 


0.07 


0.21 


-1.10 


-1.10 


-0.98 


631.303 


1.58 


0.11 


+0.27 




+0.18 


644.572 


1.00 


0.094 


-0.83 


-0.83 




699.084 


0.62 


0.050 


-1.22 




-1.44 


709.776 


0.69 


0.21 


-0.57 






710.289 


0.65 


0.065 


-0.84 




-1.06 


743.989 


0.54 






-1.18 




755.149 


1.58 






-1.36 




755.473 


0.51 






-2.28 




755.841 


1.54 






-1.47 




756.213 


0.62 






-2.71 




781.935 


1.82 


0.065 


-0.38 


-0.39 




782.292 


1.75 






-1.14 




784.938 


0.69 


0.10 


-1.30 


-1.30 




856.859 


0.73 






-2.80 




857.1085 


1.53 






-2.07 




858.421 


1.86 






-1.32 




858.787 


1.48 






-2.12 




874.958 


0.60 






-2.79 




350.1133 


0.07 


0.7 


-0.93 






357.5765 


0.07 


3.3 


-0.03 






366.3698 


0.15 


2.5 


+0.01 






588.5629 


0.07 


0.050 


-2.12 




-1.82 



B: lBiemontetalJ ((l98ll) 



VW: 
Bog: 



^nture & Wallerstein ( 2002 
Bogdanovich et al.i(il996 



A 


Elow 


EW [pm] 




log 9f 




nm 


eV 


B 


L 


B 


Bog 


L 


343.2415 


0.93 


2.1 




-0.75 


-0.51 


-0.72 


345.4572 


0.93 


1.00 




-1.34 




-1.33 


345 8940 


0.96 


1.6 




-0.52 




-0.48 


347.9017 


0.53 


2.8 




-0.69 


-1.12 


-0.67 


347.9393 


0.71 


5.1 




+0.17 


+0.12 


+0.18 


349.9571 


0.41 


2.4 




-0.81 


-1.08 


-1.06 


350.5666 


0.16 


5.1 




-0.36 


-0.62 


-0.39 


354.9508 


1.24 


1.6 




-0.40 


-0.68 


-0.72 


355 1951 


0.09 


5.8 




-0.31 




-0.36 


358.8325 


0.41 


2.7 




-1.13 


-1.25 


-1.13 


360.7369 


1.24 


1.3 




-0.64 


-0.48 


-0.70 




0.71 


3.2 




-0.60 


-0.56 


-0.58 


ItllAlll 


0.53 


3.0 




-0.93 




-0.96 


379.6496 


1.01 


1.5 




-0.83 


-1.17 


-0.89 


383.6769 


0.56 


4.7 




-0.06 


-0.22 


-0.12 


403.4091 


0.80 


0.65 




-1.55 




-1.51 


405.0320 


0.71 


2.38 


2.20 


-1.00 


-0.60 


-1.06 


408.5719 


0.93 


0.54 




-1.61 


-1.54 


-1.84 


420.8980 


0.71 


4.3 


4.26 


-0.46 




-0.51 


425.8041 


0.56 


2.6 


2.34 


-1.13 




-1.20 


431.7321 


0.71 


1.20 




-1.38 




-1.45 


444.2992 


1.49 


2.0 


2.04 


-0.33 




-0.42 


449.6962 


0.71 


3.6 


3.15 


-0.81 


-0.87 


-0.89 


511.2270 


1.66 


0.83 


0.78 


-0.59 


-0.58 


-0.85 


363.0027 


0.36 


3.8 




-1.11 




-1.11 


407.1093 


1.00 


2.1 




-1.60 




-1.66 


414.9202 


0.80 


7.5 




-0.03 




-0.04 


416.1208 


0.71 


5.8 




-0.72 




-0.59 


426.4925 


1.66 


1.50 




-1.41 




-1.63 


444.5849 


1.66 


0.83 




-1.35 






461.3921 


0.97 


2.91 




-1.52 




-1.54 


402.4435 


0.999 




1.20 






-1.13 



B:|Biemont etal.'/l98l') 

Bog: Bogdanovich et alj ( Il996h 
L: [Liung et alj ( 120061) 



(119811) . the Holweger-Miiller model with microturbu- 
lence of l.Skms"^, and the solar abundances from 
19891) . They remark that in their 



Anders & Grevesse 



metal poor stars Zr abundances from Zr I are, on aver- 
age, lower tha n those from Zr II lin es, a result similar to 
that found by iBrown et al.l (Il983h for Pop. 1 stars. For 
the analysis of the metal poor stars they used 7 Zrll 
lines and 3 of them (407. 1 414.9, and 416.1 nm ) are 
among those discarded by iBiemont et al. (1981), sug- 
gesting that these are no longer significantly blended in 
metal poor stars. They obtain for the solar Zr abundance 
2.59±0.04 from Zrl (34 lines) (7=0.22 dex, 2.53±0.03 
from Zr II (24 lines ) (7=0.14 dex, respectively. 
Bogdanovich et al. |([l£96) computed theoretical log gf. 



and used the Holweger-Miiller model with £=0.8, 
and a model computed with a code by Kipper et al.l 
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(119811) which is a modifi ed version of the AT- 
LAS 5 code ( Kurucz 1970| ). They used the damp- 
ing constants modifie d by Gald ikas ( 19881). and the 
EW measurements of iBiemont et al.l (il98lh . The re- 
sults are: A(Zr)=2.60±0.07 from Zrl (21 lines) and 
A(Zr)=2.61±0.11 from Zrll (15 lines). They adopted 
A(Zr)= 2 .60 ± 0.06. 



4. iLiung et al.l (12006.) considered only Zr II lines, they used 



log gf from experimental branching ratios measured 
by them, and radiative lifetimes taken from different 
papers. They used ID Holweger-Miiller and MARCS 
models, and a 3D model computed w ith the Stein- 
Nordlund code /S tein & Nordlundll998h for line profile 
fitting. EWs of the best 3D fitting profiles are used to 
estimate the abund ance. From 7 Z r II fine s, six of which 
in common with Biemon t et al. 1 (Il98lh . they obtain 
A(Zr)=2.58 ± 0.02 with the 3D model, A(Zr)=2.56 ± 
0.02 with MARCS model, and A(Zr)=2.63 ± 0.02 with 
the HM model. Their adopted value is A(Zr)=2.58 ± 
0.02. 

5. IVelichko et all (I2OI0I) performed the first NLTE anal- 
ysis of Zr in the solar photosphere. They inves- 
tigated all the Zr lines ana lysed by iBiemont et al 



(Il98ll) and lLjungetaP (l2006l) and selected a subsam- 
ple of two Zrl lines (424.1 and 468.7 nm) and ten 
Zrll lines (347.9, 350.5, 355.1, 405.0, 420.8 425.8 
444.2, 449.6, and 511.2nm). Analysing the iKurucz 
(l2005ah solar spectrum with a MAFAGS model at- 
mosphere (5780K/4.44/0.0) and a microturbulence of 
0.9km s^^they derive a LTE abundance of 2.33 and 2.61 
from Zr I and Zr II lines, respectively. They constructed 
a Zr model atom and derived the NLTE abundance with 
different assumptions about the rates of collisions with 
H-atoms. In their scenario, the NLTE abundance is al- 
ways higher than the LTE abundance varying in a range 
from 0.21 to 0.32 dex for Zr I and from 0.01 to 0.08 dex 
for Zrll, depending on the adopted cross sections for 
collision with H-atoms. Their best estimate of the NLTE 
abundance (correction) from Zrl and Zrll lines is 2.62 
(0.29) and 2.64 (0.03) dex, respectively. Their adopted 
Zr abundance is 2.63 ± 0.07. 



The abundance of Zr in the solar system is derived from 
the analysis of meteorites. Zirconium is not a volatile ele- 
ment, and its abundance in meteorites should agree with the 
one in the solar photosphere. To compare the abundances 
derived from the meteorites, on the cosmochemical abun- 
dance scale relative to 10® Si atoms, with the ones derived 
form the solar photosphere, on the astronomical scal e of 
10^^ H -ato ms, a couphng facto r must be derived. Lod dersI 
(l2003h and lLodders et al.l (12009 ) to Hnk the two scales in- 



troduce an average factor by looking at a sample of refrac- 
tory elements well-determined in the photosphere and at the 
same elements in meteorites. In this way the factor is not 
affected by the individual uncertainty and, if the abundance 
of one element derived from the photosphere changes, the 
conversion factor is not much affected. 



The meteor itic Zr abundance 
(l2009h . 



Lodders et al. 



is A(Zr)=2.57 ± 0.04, 
Other Zr meteoritic values that 
can be found in the literature are: 2.61 ± 0.03 accord- 

2.61 ± 0.02 



ing to 'Anders & Grevesse (1989), 2.61 ± 0.02 accord- 
ing to .Grevesse & Sauval (1998), 2.60 ± 0.02 accord- 



in g to iLoddersI (l2003h. The value of 2.53 ± 0.04 given 



in Asplund et al.^ (2009) is based on the same data from 
Lodders et al. (2009). but they give a lower value for Zr 
on the astronomical scale because they use a lower coupUng 
factor for the cosmochemical and astr onomical scales . In the 
sam e way the value of 2.5 7 ± 0.02 in lAsplund et alj (l2005h 
and Grevesse et al. (2007) is based on data from LoddersI 
(I2OO3I) . 3ut their value is lower because again a different 
scaling factor (w hich is . 03 log; units smaller than the one 
recommended in lLoddersI (l2003h ^ was used. 



4 Model atmospheres 

We base our abundance analysis on the same model at 
mospheres that we used in the previous analys is of so- 



lar abundance determinations, summarised in Caffau et al. 



(l201fll) . We rely on a time-dependent, 3D hydrodynamical 
model a tmosphere, computed with the CO^BOLD code; see 
Freytag . Steffen. & Dorch (200j) and iFreytag et alJ (I2OI0I) 
for details. This model has a box size of 5.6 x 5.6 x 
2.27 Mm^ , resolved by 1 40 x 1 40 x 1 50 grid points . Its range 
in Rosseland optical depth is —6.7 < logTRoss < 5.5. The 
complete time series is formed by 90 snapshots covering 
1 .2 h of solar time. For the spectral synthesis computations 
we selected 19 representative snapshots out of the complete 
series of 90 snapshots. 

To compute 3D-coiTections for the zirconium abun- 
dance, we make use of two ID reference atmospheres: 
the ID model obtained by horizontally averaging each 
3D snapshot over surfaces of equal (Rosseland) op- 
tical depth (henceforth (3D) model), and the hydro- 
static ID mixing-length model computed with the LHD 
code, that employs the same micro-physics and radia- 
tive transfer scheme as the CO^BOLD code, (hence- 
forth IDt.hd mod el). We define 3D-corrections as in 



Caffau et al.l(l2010l) . Agran = A(Zr)3j) -A(Zr)^3j3^, to iso- 



late the effects of horizontal fluctuations (granulation ef- 
fects), and Alhd = A(Zr)3j-, — A(Zr)j jjp, to r neasur e the 
total 3D effect. For details see ICaffau et al.1 ilOld and 
Caffau & Lu dwig (2003). 

We also consid ered the semi-empirical Holweger- 
Miille r solar model ( HolwegeJ ll967( iHolweger & Miiller 
1974 , hereafter HM) for comparison. 



5 Observed spectra 

We analysed the same four solar spectral atlases, publicly 
available (two disc-centre and two disc-integrated atlases), 
that we us ed in all our previou s works on solar abundances, 
such as in ICaffauet al.l(l201fll) . These high resolution, high 



signal-to-noise ratio (S/N) spectra, are: 
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1. the disc-integrated spectrum o 

f iKurud (12005a'), based 
on fifty solar FTS scans, observed by J. Brault and L. 
Testerman at Ki tt Peak between 1981 and 1984; 

2. the FTS atlas of' Neckel & Labsl (Il984 . observed at Kitt 
Peak in the 1980ies, providing both disc-centre and disc- 
integrated data; 

3. the disc-centre atlas of iDelbouille et alj (Il973h . ob- 



served from the Jungfraujoch. 

6 Abundance determination methods 

In principle, the best way to derive the abundance of an ele- 
ment from a spectral line is by line profile fitting. Once the 
atomic data are known, this technique allows to take into 
account the strength of the line and its shape at the same 
time. The limiting factor is that the observed solar spec- 
tra are of such a good quality that the synthetic line profile 
is often not realistic enough to reproduce the shape of the 
line in the observed spectrum. Synthetic line profiles de- 
rived from ID model atmospheres do not take into account 
the line asymmetry induced by convection, they are sym- 
met ric. In many cases NL TE effects modify the line profile 
(see lAsplund et al.ll2004l) . In case the line is not clean, but 
some contaminating lines are present in the range, the line 
profile fitting should take into account these blending com- 
ponents. Good atomic data for these blending lines are then 
necessary, and, as for the line of interest, granulation and 
NLTE effects can be limiting factors. 

To fit a line blended with other components, one has 
to optimise the agreement between synthetic and observed 
spectra also for the blending lines, allowing, in the line pro- 
file fitting process, the abundance of the blending lines to 
change as well. Even though the 3D spectrum synthesis is 
still computationally demanding, we have attempted to de- 
rive A(Zr) by properly taking into account also the blend- 
ing lines whenever possible. In the fitting procedure, the Zr 
abundance and the abundances of the blending lines are ad- 
justed independently until the best match of the observed 
spectrum is achieved (in the following referred to as method 
'A'). 

We also tested a simplified fitting approach (in the fol- 
lowing method 'B'): we synthesised one typical 3D profile 
of each blending components, in addition to a grid with dif- 
ferent abundances for the Zr line of interest. Then, in the 
fitting procedure, a profile is interpolated in the grid for Zr, 
while the profiles (more precisely the relative line depres- 
sions) of the blending lines are simply scaled, and then the 
blending lines and the Zr line are added linearly to obtain the 
composite line profile. In the fitting procedure, the Zr abun- 
dance and the scaling factors (and shifts) of the blending 
lines are adjusted until the best agreement with the observed 
spectrum is found. An advantage of this fitting procedure is 
that it works even if the atomic data of the blending compo- 
nents are poorly known. Note, however, that this procedure 
is theoretically justified only if all the lines are on the linear 
part of the curve of growth. If the lines are partly saturated. 



the line strengths of all components are underestimated, and 
so is the derived Zr abundance. Nevertheless, method 'B' 
was also applied to partly saturated lines in order to under- 
stand its limitations. 

Methods 'A' and 'B' provide the Zr abundance without 
the need to measure an equivalent width. However, we can 
formally derive the equivalent width of the Zr component 
from the synthetic line profile computed for Zr only (ignor- 
ing all blends) with the Zr abundance that provides the best 
fit of the observed spectrum. For method 'B', this is sim- 
ply the equivalent width of the Zr line that gives the best 
agreement between added-synthetic and observed profile. 

A third alternative (method 'C') is the measurement of 
the equivalent width of the line of interest, which is then 
translated into an abundance via the 3D curve-of-growth. 
For this purpose, we employ the IRAF task splotO for fit- 
ting Gaussian or Voigt profiles to the observed line. Even 
in the case that the line of interest is contaminated by 
unidentified components, which cannot be taken into ac- 
count with methods 'A' and 'B', it may still be possible 
to use splot with the deblending option to determine the 
equivalent widths of all components. This procedure may 
also be useful if the number of blending lines is too large to 
be treated conveniently by the methods ' A' or 'B' based on 
synthetic spectra. However, as method 'B', this approach is 
only justified if all involved lines are weak. 



7 Zirconium abundance from Zr i lines 



We lo oked at the Zr I lines in the sample of iBiemont et al. 
(Il98ll) . All these lines are rather weak and blended, or in 
a crowded spectral region. After inspection, we analysed 
a subsample of 1 1 lines, discarding the lines we think to 
be too heavily blended and/or in a too crowded range. We 
found our EW measurements to be affected by rather large 
uncertainties, but for the majority of the 11 lines we agree 
with the results of iBiemont et alj (Il98lb . We think that the 
cleanest fines are the ones at 480.9, 481.50, 614.0, and 
644.5 nm, but also these lines are too poor to give precise 
information on the solar abundance. We found it difficult 
to use the flux spectra to measure the EWs of so weak and 
blended lines, and decided to use only the disc-centre spec- 
tra, which are sharper From the four clean lines mentioned 
above, we obtain, A(Zr)= 2.65 ± 0.07. If we remove also 
the 644.5 pm line we obtain A(Zr)= 2.62 ± 0.04. This re- 
sult is in perfect agreement with the abundance obtained 
from Zrll lines (see below), but we consider this agree- 
ment fortuitous. From the IDlhd model we find a slightly 
highervalueof A(Zr)=2.67±0.07, andfromthe (3D) model 
A(Zr)=2.74±0.08. 



■ |http : //iraf ■ noao ■ edu/| 



www.an-joumal.org 



© 2010 WILEY-VCH Verlag GmbH&Co. KGaA, Weinheim 



794 



E. Caffau el at.: Zirconium abundance in the Solar photosphere 



<3D> 
1D, „n 



^1.4 
1.2 

-i.or 

-0.8S, 
; o 

^0.6| 

LU 

-0.4-0 

-0.2 
:0.0 



-2 

log T 



1 



Fig. 1 Temporal and horizontal average of the tempera- 
ture profile of the 3D model (solid line) and the temperature 
profile of IDlhd model, shown as a function of (Rosseland) 
optical depth. In addition, the equivalent width contribution 
functions (lower, roughly Gaussian-shaped curves) for the 
405.0 nm Zrll line are plotted on the same optical depth 
scale for disc-centre (solid line) and disc-integrated spectra 
(dashed hne). 

8 Investigation of Zr ii lines 



We picked up all the Zr II lines of lBiemont et a l. Vl98l') and 
the line at 402.4pm analysed in iLiung et al.l C2006.) . Our 
starting sample was formed by 32 Zr II lines. After inspec- 
tion of the observed data, and comparison with ID synthetic 
spectra, we discarded 1 1 lines that we think are too heavily 
blen ded. After analysis, \y e discarded four other lines that 
also iBiemont et al.l (Il98ll) rejected because they provided 
a too high ab undance. Ou r final sample is formed by the 
seven lines of iLiung et al.l (120061) . that we discuss one by 
one in the next section, and 10 lines that are used for the 



abundance analysis in iBiemont et al.l(ll981h . 

These Zr II fines have lower level energies below 1 .7 eV. 
If we look at the equivalent width contribution functions for 
disc-centre of the 3D model, we can see that they are formed 
in a range in tKoss between —4 and 0, with the maximum 
in the range —1.0 < tross < —0.5 (see Fig.[T]i. Only two 
lines (350.5 and 355.1 nm) are formed over a larger range, 
with a maximum around tross ~ —2.0, because these two 
lines have smaller lower level energies, and are significantly 
stronger than the others. 

The abundance results we obtained from the Zr II lines 
are summarised in Table|3] We averaged the EW measured 
in the two disc-centre and the two disc-integrated spectra, 
respectively. The last column of the table indicates if the re- 
sult refers to disc-centre or disc-integrated spectra. From the 
complete sample of 17 lines we obtain A(Zr)=2.612±0.093, 
with 2.638 ± 0.078 and 2.586 ± 0.100 from disc-centre and 
disc-integrated spectra, respectively. The line at 347.9 nm 
gives a very low abundance, while the line at 355. Inm 
gives a rather high one. Both these two lines lie in a very 
crowded region of the spectrum and it is very difficult to 



place the continuum. When we remove these two lines 
from the sample, we have 15 Zrll lines and we obtain 
A(Zr)=2.616±0.060, with 2.643 ± 0.050 and 2.589 ± 0.057 
from disc-centre and disc-integrated spectra, respectively. 
We note that the 3D abundance from flux is systematically 
lower than the one from the intensity spectra. 



W hen we look only at the seven lines from lLjung et al. 
('2006'), we find that the line-to-line scatter is reduced with 
respect to the sample of 15 lines. The 3D abundance we 
find is 2.653 ± 0.022 and 2.623 ± 0.031 for disc-centre and 
disc-integrated data, respectively. The abundance averaged 
over both disc-centre and disc-integrated spectra is 2.638 ± 
0.031. We find a very low line-to-line scatter, but again the 
disc-centre spectra give a systematically higher abundance 
than the disc-integrated spectra, by an amount comparable 
to the line-to-line scatter of w 0.03 dex. 

We found similar results for Fe, Hf, and Th, where the 
abundances derived from disc-integrated spectra are lower 
by 0.02, 0.05, and 0. 02 dex, respectively , th an the disc 



centre abundances (see ICaffau et al.ll2008l and ICaffau et al. 



1201 0" for details). A similar behaviour is seen when using 
ID model atmospheres. We have presently no obvious ex- 
planation for this result. 

The zirconium abundance derived from ID models de- 
pends on the choice of the microturbulence parameter ^. The 
microturbulence we find for the sample of 15 Zrll lines, 
by requiring that the relation between EW and A(Zr) ob- 
tained with the ID model has the same slope as that of the 
3D model (method 3a in Steffe n et al. ( 2009 )) , is c onsis- 
tent with the results obtained bv lSteffen et al I ( 2009b from 



a sample of Fell lines. We adopt f = 1.0 km s ^ for the 
analysis of disc-integrated spectra, and obtain A(Zr)^3j3) = 
2.597 ± 0.053 and A(Zr)LHD = 2.575 ± 0.057. For disc- 
centre data, we choose ^ = 0.7 km s^^ and the result is 
A(Zr)(3D) = 2.652 ±0.044 and A(Zr)LHD = 2.635±0.045 
(for ( = l.Okms-i we get A(Zr)(3D) = 2.602 ± 0.067, 
A(Zr)LHD = 2.587 ± 0.072). 

3D corrections are small in absolute value and depend 
on the value of the microturbulence. The average 3D cor- 
rection with respect to the reference IDlhd model is pos- 
itive, Alhd = -1-0.011 ± 0.012 dex, while, with respect to 
the (3D) model it is slightly negative, Agian ~ —0.008 ± 
0.016dex, again taking ^ = 0.7kms^^ for disc-centre and 
^ ~ 1.0 km s^^ for disc-integrated spectra. 

With the HM model we obtain A(Zr)HM = 2.679 ± 
0.044 and A(Zr)HM = 2.619 ± 0.053 for disc-centre and 
integrated disc, respectively. Taking into account both lines 
of disc-centre and integrated disc, we obtain A(Zr)HM = 
2.649 ± 0.057, and when applying the Agran correction, we 
have A(Zr)HM = 2.641 ± 0.057. 

For our final determination of the solar zirconium abun- 
dance, we adopt the complete sample of 15 Zrll lines, and 
our recommended value is A(Zr)3D = 2. 62±0. 06, obtained 
with 3D model atmosphere. 
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Table 3 The results from the 7 lines of Zr II from lLiung et al. I (l2006h . and other 10 lines from lBiemont et alj (Il981h 
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9 Analysis of a selected subsample of Zr ii 
lines 

We particu l arly i nvestigated the Zrll lines selected by 
Ljung et al.l (120061) for the solar abundance determination. 
We find this sample interesting because it contains: two 
mostly clean lines (A 405.0, 420.8 nm), a line blended on 
the blue side (A 449.6 nm), two lines blended on the red 
side (A 402.4, 425.8 nm), and two lines blended on both 
sides (A 444.2, 5 11. 2 nm) 

9.1 The Zr ii line at 402.4 nm 

This weak Zr II line is on the blue wing of a strong blend of 
Ti I and Fe II where the Ti I line is the dominant component. 
On the blue wing of this Ti-Fe blend, there is also a Ce II 
line. Close by, on the red side of the Ti I line, there is a 
strong line of iron blended with a much weaker Nd II line. 



The zirconium abundance is determined with method 'B' 
(see Sect.|6]i: to take into account the effects of the lines 
present in the range of the Zrll line, we computed a 3D 
profile for each of the Ce II line at 402.44, the Ti I line at 
402.45 nm, and the strong Fe I line at 402.47 nm, to be used 
in the fitting process. The comparis on of the best fitt i ng 3D 
synthetic profile and the observed iDelbouille et"al] d 1973b 
disc-centre spectrum is shown in Fig.|2] The Zr abundances 
and EWs resulting from the best fit of the four observed 
spectra are listed in Table[3] 

The situation of the 402.4 nm line may be compared 
to that of the 425.8 nm Zrll fine (see Sect.|94]l, and the 
444.2 nm Zr II line (see Sect l9.5b . for which we show below 
that the abundance is underestimated by about 0.03 dex and 
0.01 dex, respectively, when using the approximate method 
'B'. 

The use of splot for the EW measurements is difficult 
in this case, also with the deblending option. In fact, several 
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402.43 402.44 402.45 402.46 402.47 402.48 
X[nm] 

Fig. 2 The observe d solar disc-centre spectrum of 
Delbouille et al. (1197 3') in the region of the Zrll line 
at 402.4 nm (black symbols), superimposed on the best 
fitting 3D synthetic profile obtained with method 'B' 
(solid green/grey). The graph near the bottom shows the 
synthetic-observed residuals. The Zr abundance derived 
from this fit is A(Zr)=2.602. 



lines should be introduced in the range, making the profile 
fitting unstable. 

We also fitted the line profile with a grid of 3D syn- 
thetic spectra computed with zirconium only. We limited 
the fitting to the blue wing and the line core, neglecting 
the red wing, which is blended with the other lines (see 
Fig EJ. In this kind of fit, neglecting the contribution of the 
other lines, the continuum placement, when considered as a 
free parameter, is problematic. When fixing the continuum, 
the abundance from the fit is in close agreement with the 
abundance determination from method 'B', in the case of 
the disc-centre spectra within 0.03 dex (cf. Fig.[3]l. For the 
disc-integrated spectra we find the fit to be ambiguous. The 
abundance comes out larger than from method 'B' by about 
0. 1 de x, but the fit is convin cing only when taking the at- 
las of Neckel & Labs ( 19841) . although even in this case the 
core of the line is not well fitted. We find a different result 
from the two disc-integrated and the two disc-centre spec- 
tra, respectively, that we cannot easily explain. Maybe the 
differences are related to the fact that lines in flux spectra 
are more broadened than at disc-centre. 

9.2 The Zr ii line at 405.0 nm 

This line seems to be unblended. A weaker line is on the red 
side at about A A 0.016 nm from the centre of the Zr II line, 
and a stronger Fe I line is also on the red side at about AA 
0.035 nm. Looking at the spectra one could conclude that 
the line is clean and that the abundance can be deduced from 
line profile fitting with a synthetic grid computed with only 
Zr and/or by easily measuring the EW. When we fit the line 
profile with a grid that takes into account only zirconium, 
the result is very satisfactory for all observe d spectra, except 
perhaps for the disc-integrated spectrum of Neckel & Labsl 
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Fig. 3 In the top panel, the obse rved solar disc-centre 
spectrum of iDelbouille et al.l (Il973h in the region of the 
Zrll line at 402.4 pm (black symbols), superimposed on 
the best fitting 3D synthetic profile (solid green/grey) ob- 
tained with a grid computed with only zirco nium. In this 
fit A(Zr)=2.644. From the fit of the Neckel & Labsl (Il984 
disc-centre profile, we get A(Zr)=2.669. The lower panel 
shows a zoom in the fitting range. 



(119841) (see 'NF', Fig.|4]l. The results are in close agreement 
with the abundance determinations from the EWs measured 
with splot. We also compared the theoretical line bisec- 
tor from the 3D sy nthetic line profi l e to th e observed disc- 
centre spectrum of IDelbouille et al.l (Il973h . and we find an 
excellent agreement (see Fig.|5l). The shift in wavelength be- 
tween observed and syntheti c line bisector is mainly due to 
the fact that the spectrum of IDelbouille et aL ( 1973 ) is not 
corrected for gravitational redshift. 

9.3 The Zr ii Une at 420.8 nm 

This line is strong and sensitive to the adopted damping 
constants. The red side is more or less clean, while on the 
blue side there is a weak line, and close by another weak 
line. The line profile fitting is done with method 'A', us- 
ing Zr only, and restricting the fit to the wavelength range 
420.891 . . . 420.913 nm in order to avoid the weak blends 
in the far blue wing of this Zr II line. The continuum level 
is fixed in the fitting proced ure. The fittin g of the disc- 
integrated line profile in the iKuruczl ( 2005a ) spectral atlas 
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Fig. 4 The observed solar disc-integrated (top) and disc- 
centre (bottom) spectra in the region of the Zrll Une at 
405.0pm (black symbols), superimposed on the best fitting 
3D synthetic profile (solid green/grey) obtained with a grid 
computed with only zirconium. 
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Fig. 5 The observed solar spectrum of iDelbouille et al. 



(Il973h for the the Zrll line at 405.0pm (dashed black), to- 
gether with its line bisector (plotted on a ten-times expanded 
abscissa), superimposed on the 3D synthetic profile (solid 
black) with its line bisector. 
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Fig. 6 The o bserved solar disc-integrated spectrum of 



Kuruca (12005 al) (black symbols) in the region of the Zr II 
line at 420.8 pm, superimposed on the best fitting profile 
(solid green/grey) derived from a grid of 3D synthetic spec- 
tra computed with only zirconium. The fit corresponds to 
A(Zr)=2.58. 



9.4 Tlie Zr ii line at 425.8 nm 

This Zr ll line is on the blue wing of a much stronger Fe II 
line that has close by, on the red side, a very strong Fel 
line, blended with a weaker Fe II line (see Fig.|7ll. We tried 
different fitting procedures. 

First, we fitted the Zr II line profile with a grid of 3D 
synthetic profiles computed with only Zr, restricting the fit- 
ting to t he blue wing of th e Zrline. For the disc-centre spec- 
ti-um of IDelbouille et al.1 (Il973 ^ we obtain A(Zr)=2.65 ± 
0.01. The fit looks quite satisfactory (see Fig.|7l), but it is 
not clear how much the strong lines close by affect the Zr II 
line, and hence by how much A(Zr) is overestimated. 

Second, we applied method 'B', taking into account 
the contribution of the neighbouring iron lines as a lin- 
ear combination of scaled 3D synthetic profiles, and obtain 
A(Zr)=2.61 ± 0.01. 

Finally, we also fitted the complete range with method 
'A', properly taking into account the three iron lines, with 
A(Zr) and A(Fe) as free fitting parameters. We obtain a 
Zr abundance of 2.643 ± 0.001 and 2.627 ± 0.000 from 
disc-centre and disc-integrated spectra, respectively. Even 
though the iron lines are not very well reproduced (see 
Fig.[8]l, the Zrll line is fitted very well and is not affected 
by the minor deficiencies of the iron lines. The equivalent 
widths given in Tableware derived from the abundances ob- 
tained with this fitting method. 



is shown in Fig.|6] The abundance derived from line pro- 
file fitting and from the EW measured with splot are in 
very good agreement (within 0.0 07 d ex) for t he obse rved 
spectra o f Delbouille et al.l (119731) and lKurucj (l2005al) . For 
Neckel & Labsl (Il984j) . the abundance from line profile fit- 
ting is lower with respect to the abundance derived from 
the EW measurements by 0.04 dex and and 0.01 dex, for the 
disc-integrated and disc-centre spectrum, respectively. 



9.5 Tlie Zr ii line at 444.2 nm 

This line is surrounded by two stronger features, an iron line 
on the blue side at 444.283 1 nm and a blend of three lines of 
iron on the red side. Of these three Fe I lines on the red side, 
the one at 444. 3194 nm is much stronger than the other two 
components. 

Clearly, the Zr II line is affected by the two Fe I features 
on both sides, and cannot be treated in isolation. Since the 
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Fig. 7 The observe d solar disc-centre spectrum of 



Delbouille et alJ (119731) in the region of the Zrll line at 
425.8pm (black symbols) is superimposed to the 3D best 
fit (solid green/grey), obtained with a grid of 3D synthetic 
profiles computed with only Zr. 
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Fig. 8 The observed iNeckel & Labs (11984 spectrum in 
the region of the Zr II line at 425.8 pm (black symbols) is su- 
perimposed to the best fitting 3D synthetic spectrum (solid 
green/grey) obtained with method 'A'. Upper panel: disc- 
integrated spectrum, lower panel: disc-centre spectrum. 

two iron features strongly influence the shape of the Zrll 
line, the Zr abundance cannot be derived from this line by 
fitting with a grid of 3D synthetic profiles computed with 
zirconium only. 

The most reliable Zr abundance is obtained from the full 
fitting, method 'A': we fitted the observed Fe-Zr-Fe blend 
with A(Zr) and A(Fe) as free parameters. The Zr abundance 
derived in this way is 2.673 ± 0.004 and 2.650 ± 0.005 for 
the disc-centre and disc-integrated spectra, respectively (cf. 
Fig.|9|l. The equivalent widths given in Table[3]are derived 
from the abundances of this fitting method. 

For completeness, we also derived the Zr abundance 
with methods 'B' and 'C. 

For method 'B', the profiles of the four iron blending 
features are represented by scaled 3D synthetic Fe line pro- 
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Fig. 9 The observ ed solar disc-centre spectrum of 



iNeckel & Labsl (Il984 in the region of the Zrll line at 
444.2 pm (black symbols) with the 3D synthetic best- 
fit-profile obtained with method 'A' superimposed (solid 
green/grey). Upper panel: disc-integrated spectrum, lower 
panel: disc-centre spectrum. The residual is below each plot. 

files, which are then added to the profile of the Zr II line in 
order to obtain an approximate profile of the complete blend 
(ignoring saturation effects). With this approach, the best 
fits are obtained with Zr abundances of A(Zr)=2.67 ± 0.01 
and A(Zr)=2.64 ± 0.01 for disc-centre and disc-integrated 
case, respectively. The fit of the observed disc -integrated 
spectrum is somewhat better that in the case of the disc- 
centre observations. 

In method 'C, A(Zr) is found from the equivalent width 
of the Zrll line, as obtained from splot with the deblend- 
ing option. Representing the Fe feature on the red side of the 
Zr II line either by a single or by three Voigt profiles turned 
out to have no influence on the resulting zirconium abun- 
dance. We find A(Zr) derived in this way to be very close 
to the result of method 'B' in the case of the disc-centre ob- 
served spectra. In the disc-integrated case, A(Zr) is about 
0.057dex smaller than what is obtained with method 'B'. 

9.6 Tlie Zr ii line at 449.6 nm 

This Zrll line (EW w 3pm) is contaminated by a much 
stronger Cr I line (EW « 9 pm). Certainly, saturation effects 
cannot be ignored in this case. 

Based on the full 3D line profile fitting of the complete 
blend using method 'A', with A(Cr) and A(Zr) as free fitting 
parameters, we obtain a Zr abundance of A(r)=2.660±0.003 
for disc-centi-e (cf. Fig.IIB, and A(Zr)=2. 642 ±0.005 for in- 
tegrated disc. These Zr abundances correspond to equivalent 
widths of 3.1 1 and 3.41 pm respectively (see Table[3]l. 

The fit resulting from method 'B', scaling the Cr pro- 
file and interpolating in the grid of 3D synthetic Zr pro- 
files to achieve the best agreement between synthetic and 
observed profile, is shown in Fig.[TO] for the disc-centre 
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Fig. 10 The observ ed solar disc-centre spectrum of 
Neckel & Labsl (Il984 in the region of the Zrll Hne at 
449.6 nm (black symbols) with the best fitting 3D synthetic 
profile obtained from method 'B' with A(Zr)=2.63 superim- 
posed (solid green/grey). The graph in the lower part shows 
the synthetic — observed residuals. 



case. It is of similar quality as the one obtained from the 
full fitting procedure with method 'A' (Fig.[TT)- Method 'B' 
gives A(Zr)=2.63 and 2.60 for disc-centre and integrated 
disc, respectively. As expected, the approximate method 'B' 
tends to underestimate the Zr abundance, in this example 
by 0.027 and 0.042 dex for disc-centre and integrated disc, 
respectively. The average equivalent width we derive with 
this method for the observed disc-centre and disc-integrated 
spectra are EW=3.01 and 3.25pm, respectively, 3.2 and 
4.7% lower compared to method 'A'. 

When we apply method 'C using splot with the de- 
blending option, we cannot reproduce the Cr line profile 
with a single Voigt profile, but we can formally have a good 
reproduction of the profile with two Voigt profiles, but this 
option cannot take into account properly the wings of the 
Cr II hne. 



9.7 The Zr ii line at 511.2 nm 

For this line, the placement of the continuum is difficult, 
but the uncertainty induces an error in the EW measure- 
ment of less than 4% and hence a difference in the Zr abun- 
dance of at most 0.02 dex. The line shape suggests that it 
is blended on both wings, and in fact three weak blend 
lines are expected according to the Kurucz line listd; two 
Kl lines at 511.2212 and 5 11. 2246 nm, and a Smll line at 
511. 2294 nm. In addition, several weak molecular lines are 
expected in the range (see Fig.[T2]i. 

The average EWs derived with method 'B' for disc- 
centre and integrated disc, respectively, are 8.4 and 9.2pm, 
corresponding to zirconium abundances of A(Zr)=2.69 and 
2.67. When measuring the EW with splot (method 'C'), 
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Fig. 11 The observed solar disc -centre spectrum of 



Neckel & Labsl (Il984h in the region of the Zrll hne at 
449.6 nm (black symbols) with the best fitting 3D syn- 
thetic profile obtained from method 'A' superimposed (solid 
green/grey). The graph in the lower part shows the synthetic 
— observed residuals. 
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Fig. 12 The observ ed solar disc-centre spectrum of 



' |http : / /kurucz ■ harvard ■ edu/linelists ■ html| 



Delbouille et al.l (Il973h in the region of the Zrll line at 
511.2nm (black solid) with the identifications of the lines 
in the range according the line list of Kurucz. 

we find similar results, but the uncertainty in the measure- 
ments is larger, giving an indetermination in the abundance 
determination of about 0.06 dex. 

10 Conclusions 

We analysed the zirconium abundance in the solar photo- 
sphere, investigating a selected sample of lines of both Zr I 
and Zrll. The Zrl lines are weak and heavily blended so 
that only four of them are acceptable for abundance work. 
However, the present analysis of the zirconium abundance 
relies primarily on 15 lines of Zrll that we found suitable 
for this purpose. 

We have applied three different fitting strategies to de- 
rive the abundance from Zr lines that are blended by lines 
of other elements. In the case that all components making 
up the blend are weak, the different methods give consistent 
abundances. If, however, stronger lines are involved (for an 
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example see Sect. |9.6l l. the methods that ignore saturation 
effects may severely underestimate the Zr abundance by up 
to 0.05 dex, even though the result of the fitting may appear 
pleasing to the eye, and the reduced may be close to one. 

We find a good agreement between A(Zr) derived from 
the Zr I and the Zr II lines, but, due to the scarcity of Zr I 
lines, we consider this result as fortuitous. The abundance 
from disc-centre spectra is systematically higher than the 
one from disc-integrated observations. A similar result is 
also found for other heavy elements like Fe, Th, Hf, both 
with 3D and ID model atmospheres. This may indicate a 
problem with the thermal structure of the models, rather 
than a physical abundance gradient in the solar atmosphere. 
Further investigations are necessary to find an explanation 
for this small discrepancy. 

Our recommended solar zirconium abundance is based 
on the 3D result for 15 Zr II lines, and is A(Zr)= 2.62±0.06, 
where the uncertainty is the line-to-line scatter of the se- 
lected sample of Zrll lines. This value is at the upper 
end of the solar zirconium abundances found in the liter- 
ature. Still, within the mutual error bars, this result is in 
good agreement with the meteoritic zircon ium abundance 
of A(Zr)=2.57 ± 0.04 dLodders et al. Il2009h . 
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